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Abstract 

We establish a general gluing theorem for constant mean curvature solutions of 
the vacuum Einstein constraint equations. This allows one to take connected sums 
of solutions or to glue a handle (wormhole) onto any given solution. Away from this 
handle region, the initial data sets we produce can be made as close as desired to the 
original initial data sets. These constructions can be made either when the initial 
manifold is compact or asymptotically Euclidean or asymptotically hyperbolic, with 
suitable corresponding conditions on the extrinsic curvature. In the compact setting 
a mild nondegeneracy condition is required. In the final section of the paper, we 
list a number ways this construction may be used to produce new types of vacuum 
spacetimes. 

1 Introduction 

1.1 The constraint equations and surgery 

The Einstein equations for the gravitational field on a Lorentzian four-manifold Z form, 
modulo diffeomorphisms, a locally well-posed hyperbolic system [11], and vacuum solu- 



tions may be obtained from a set of Cauchy data on a three-dimensional spacelike slice 
E C Z . A vacuum initial data set for this problem consists of two symmetric tensors 7 
and IT on E, which correspond respectively to the induced metric and second fundamental 
form on E. Thus 7 is positive definite while II is at least apparently unrestricted. How- 
ever, in order for the Einstein evolution to exist, at least for a short time, it is necessary 



and sufficient [LI] that these tensors satisfy the constraint equations: 



divTI 



vtrn 

2 



|n|; + (trn 



(i) 

(2) 



Here V is the Levi-Civita connection for 7, i? 7 its scalar curvature, and the divergence, 



div 7 , trace, tr 7 
respect to 7. 



and norm squared of n, | H| . 



n a bH c d'y ac "f bd , are all computed with 
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Having fixed the three-manifold E, it is not apparent that solutions (7, II) of these 
constraint equations exist, or how to find them. However, there is a systematic procedure 
which generates solutions in a large number of cases, which we review now. Let us begin 
with a pair of tensors (7, n), which do not necessarily satisfy the constraint equations, 
but which have the property that the mean curvature r = tr 7 n is constant. The method 
we are describing relies crucially on this assumption. Decompose n into its trace-free 
and pure trace parts: 

n = /i + ~ 7 . (3) 

The first of the constraint equations (||) requires that div^ = 0, since r is constant, so 
that fj, is what is usually called a transverse-traceless (TT) tensor field. The idea is to 
modify the pair (7, n) by changing 7 and the trace-free part /i of n by a conformal factor; 
by a judicious choice of this factor, the new tensors will satisfy the constraint equations. 
More specifically, for any <f> > we set 

7 = </> 4 7 , (4) 
ft = /i + ^7 = fV + ^ 4 7- (5) 

The factor in fl is the only one for which fx is divergence-free with respect to 7. 
(Note that II and fl as they appear in (||) are covariant tensors; the factor i^ -10 replaces 
(f>~ 2 if we work with II and fl in contravariant form.) This new pair (7,1!) satisfies the 
constraint equations provided the conformal factor (f) satisfies the Lichnerowicz equation 

- ±R4> + \\fi\ 2 r 7 ~ ^r 2 cf> 5 = 0. (6) 

Here the Laplacian, scalar curvature and norm-squared of (i are all computed relative to 
7. By well-known analytic techniques one may find solutions of this semilinear equation, 
and so altogether this procedure produces many admissible pairs of tensors satisfying 
the constraint equations. In effect, the constancy of the mean curvature r decouples 
the equations (Q), @, reducing them to the facts that \x is transverse-traceless and 4> 
satisfies the Lichnerowicz equation. These formulations are equivalent and so we label 
initial data sets in one of three equivalent ways: as (7, n) if these tensors satisfy the 
constraint equations, or as (7, n, <f>) if (f> is the conformal factor which alters (7, n) to 
tensors satisfying the constraint equations, or finally (suppressing r) as (7,/i,^). When 
<j) is omitted in this notation, then <f> = 1 is implied. 

When X is compact, complete existence results are available in this CMC case [18|; 
analogous results are known when X is either asymptotically flat [0, [15], ||, f8[or 
asymptotically hyperboloidal (3|, In the non-CMC case, when r is no longer constant, 
only partial results have been obtained [19], [12], 14], [20] but this direction warrants 



greater attention since the existence of constant mean curvature slices is apparently a 
restrictive assumption in relativity ||, ||. 

In this paper we address the following question. Suppose we are given a three-manifold 
X along with an initial data set (7, n) on it which satisfy the constraint equations. Then 
is it possible to modify X by 'surgery' to obtain a new, topologically distinct manifold X', 
and to find tensors 7' and n' on X' which themselves satisfy the constraint equations? 
Since we can find solutions by applying the procedure above to X' in the first place, the 
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more precise question we propose to study is whether it is possible to perform the surgery 
on a geometrically small subset of E and to find solutions (7', II') which are very close 
to the initial solutions (7, II) away from this subset? 

Let us phrase this more concretely. There are many different types of surgeries 
possible on higher dimensional manifolds, and really only two types in three dimensions. 
One is akin to Dehn surgery, while the other, on which we focus, consists of 'adding a 
handle' in the following sense. Choose any two points Pi,p% £ E and let B\ and B2 be 
small (metric) balls centered at these points. Excise these balls to obtain a manifold 
E" = E \ (B\ U B2). The boundary of E" consists of two copies of S 2 , the same as the 
product of an interval with a two-sphere, I x S 2 . Identifying these two manifolds along 
this common boundary yields a new manifold 

s = s" y /x s 2 . 

{±1}xS 2 

Informally, we propose to restrict the original initial data set (7,1!) to E", modify this 
data slightly, and then find a one-parameter family of extensions (7 e ,II £ ) of this data 
to E, e £ (0, £o)) each of which also satisfy the constraint equations. The parameter e 
measures the size of the handle. We shall do this in such a way that on any compact 
set K away from the handle, (7 e ,II e ) converges to (7,11) as some power of e. The data 
on the handle will converge to zero. In the compact setting, the only condition needed 
to make this procedure work is a very mild nondegeneracy condition as well as the 
assumption that II ^ 0; the precise statement is contained in the main theorem below. 
In particular, we cannot handle compact 'time-symmetric slices', i.e. those with IT = 0, 
in which case the constraint equations reduce to the vanishing of scalar curvature, and 
known topological obstructions (see [ |32| for example) prevent this sort of construction 
from working in general then. 

There are two main cases of this construction. In the first, E consists of two com- 
ponents, Si and E2, and the two points pi,P2 he in these separate components. Then 
we are producing solutions of the constraint equations on the connected sum Ei^S2 of 
these two components. In the other case, E is connected, and this construction may 
be regarded as demonstrating that a 'wormhole' can be added to most solutions of the 
constraint equations. Both cases are of clear physical interest. 

The analytic techniques used in the proof here are in the spirit of many other gluing 
theorems for a variety of other geometric structures. Gluing theorems are fundamental 
in gauge theory, and similar techniques have been developed and applied to metrics 
of constant scalar curvature [29], p(|, minimal and constant mean curvature surfaces 



[21], [p7|l , [28], and a variety of other geometric problems. Another recent result using 



the same circle of ideas, and relevant to this area of general relativity, is contained in 



[16]; this paper proves the existence of time-symmetric asymptotically flat solutions to 
the constraint equations (i.e. scalar flat asymptotically Euclidean metrics) which evolve 
to spacetimes which are identically Schwarzchild near spatial infinity. The technique 
common to many of these papers consists of a two-step procedure. First a one-parameter 
family of approximate solutions to any one of these problems is constructed. These are 
chosen so that the error measuring the discrepancy of these approximate solutions from 
exact solutions tends to zero, and so the strategy is to perturb these approximate solutions 
to exact solutions using a contraction mapping or implicit function theorem argument. 
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The complication here is that as the parameter e converges to zero, the geometry is 
degenerating and some care must be taken to control the linearization of the nonlinear 
equation associated to the geometric problem. A slightly different and ultimately simpler 
approach is contained in |27]], p§| ; however, we shall follow the earlier and more familiar 
route here because in the present context either method requires approximately the same 
amount of work. 

In the remainder of this section we set up some notation, give some definitions, state 
precise versions of our main results, and provide a more detailed guide to the rest of the 
paper. 

1.2 Statement of main results 

We fix a three- manifold E along with a conformal class [7] on £. For any two points p\,P2 
on this manifold, we refer to the triple {[y],Pi,P2) as a marked conformal structure on 
£. We have already alluded to a nondegeneracy condition required for the construction: 

Definition 1 A marked conformal structure ([y],Pi,P2) on a compact manifold £ is 
nondegenerate if either of the following situations is true: 

• E is connected and there are no nontrivial conformal Killing vector fields on £ 
which vanish at either p\ or P2, or 

• The points p\ and P2 lie on different components of £ and there are no confor- 
mal Killing fields which vanish at pj but which do not vanish identically on the 
component of £ containing pj, for j = 1,2. 



Remark 1 To avoid trivialities, we always assume that either E is connected or else has 
two components, each of which contains one of the points pj . Note that this nondegeneracy 
condition is very mild and holds for any fixed conformal class [7] when the points pj 
are chosen generically. In the asymptotically Euclidean and asymptotically hyperboloidal 
cases, which are the only noncompact cases of interest to us here, this nondegeneracy 
condition is not needed. We discuss these cases in more detail in §7. 

Assume that E is compact, and fix two points pi, P2 on it. Suppose that (7,11) is 
an initial data set with II ^ and r = tr 7 II constant, so that in particular the two 
constraint equations (g) are satisfied. It is often more convenient to regard (7,/^) as 
the initial data set instead, where [i is the transverse-traceless part of II. For any (small) 
R > 0, fix balls Bj = Bji(pj) of radius R (with respect to 7) around the points pj, and 
define 

Efc = £\(5iUB 2 ). (7) 

We also let Eq = E* = E \ {pi,P2}- We shall modify the metric 7 conformally in each 
of the punctured balls Bj \ {pj} to obtain a metric j c on £* with two asymptotically 
cylindrical ends. These ends have a natural parameter t = — logr, where r is the geodesic 
distance in Bj to pj. According to (||) the transverse-traceless tensor fi changes to a 
tensor /u c which is transverse-traceless with respect to j c . This gives an asymptotically 
cylindrical pair (7 C ,/^ C ) together with a function ip on E* which solves the associated 
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Lichnerowicz equation @ (with coefficients determined by this pair of tensors) . We also 
refer to the triple (j c , // C jV0 as an initial data set. 

The construction proceeds as follows. We first identify long pieces, of length T, of 
these two cylindrical ends with one another to obtain a family of metrics jx on the 
new manifold Ey obtained from E by adding a handle. For each T, St is of course 
diffeomorphic to the manifold E described in the previous section. This process involves 
cutoff functions and is not canonical. We also use cutoffs to patch together the values of \x c 
on these ends to obtain fix, and similarly an approximate solution ipx of the Lichnerowicz 
equation. We write this out more carefully in §|2[ Although it is not hard to construct \ix 
so that it still has vanishing trace relative to 77 1 , it is no longer necessarily true that \ix 
will be divergence-free. The triple (77 1 , \iti^t) does not solve the constraint equations, 
but instead gives an approximate solution, with an error term tending to zero as T —* 00. 
We show that it is possible, for T sufficiently large, to correct this data to obtain exact 
solutions of the constraint equations. We do this in the following order. First fix is 
altered by a small correction term ax, using a well known second order linear elliptic 
system (usually called the vector Laplacian) derived from the conformal Killing operator 
on (St, 7t), to obtain a tensor fix = fJ-x — o'x which is transverse-traceless with respect to 
7t- This involves a careful analysis of the mapping properties of this operator with good 
control as T — > 00 and is the topic of §||. We then use the conformal method, and more 
specifically the Lichnerowicz equation (||), to modify the function tpx to some ipx, so that 
(7t> MTj ipr) is an initial data set, i.e. gives an exact solution of the constraint equations. 
This involves first estimating the error terms for this approximate solution, which is 
done in §^j. We next analyze the linearization of the Lichnerowicz equation uniformly as 
T — > 00. In particular, we introduce a class of weighted Holder spaces C^' a (Y,x) which 
correspond to the standard Holder spaces away from the handle and determine the range 
of weights 5 for which the inverse of the linearized Lichnerowicz equation is uniformally 
bounded independent of T. This is carried out in §|S[ Finally, in §|6] we use a contraction 
mapping argument to produce the exact solutions by finding a function rfx in the ball 
M v of radius ve~ T l A in C^ +2,q (St) for 5 G (0,1), so that ipx = ?Pt + ffr satisfies the 
Lichnerowicz equation. This establishes our main theorem which we now state (in the 
form appropriate for initial data on closed manifolds). 

Theorem 1 Let (S, 7, H;p\ ) p2) be a compact, smooth, marked, constant mean curvature 
solution of the Einstein constraint equations; thus n = fi + 5T7 with r constant and [i 
transverse-traceless with respect to 7. We assume that this solution is nondegenerate 
and also that H 0. Then there is a geometrically natural choice of a parameter T 
and, for T sufficiently large, a one-parameter family of solutions (T,x,Tx,Hx) of the 
Einstein constraint equations with the following properties. First, the three-manifold Sy 
is constructed from S by adding a handle, or neck, connecting the two points p\ and P2 ■ 
Next, 

T T = 4> T It and U T = V^ 2 fix + -^t^xIt, 

where ipx = ipT + f]T and fix = [i — o~x are very small perturbations of the conformal 
factor ipx and approximate transverse-traceless tensor fix defined informally above and 
more carefully in §2 below. The tensor fix is transverse-traceless with respect to Tx 
and the function ipx satisfies the Lichnerowicz equation with respect to the pair (7t,/^t)- 
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Finally, the perturbation terms satisfy the following estimates as T — > oo: 

• On SJj, 

|r T - 7 | 7 < Ce" T/4 , |n T - n| 7 < CT 3 e~ 3T/2 ; 

• On i/ie nec& region, there is a natural choice of coordinate system (s,9), —T/2 < 
s < T/2 and 9 E S 2 , such that if ho is the standard metric on S 2 , then 

7T = ds 2 + h + 0(e~ T/2 coshs), 

^ T = 2e~ T/i cosh(s/2) + 0(e" T/4 ), |r? T | < Ce^ 1+S)T/4 {cosh{s/2)) 5 
for some weight 5 £ (0, 1), and finally 

|/2 T | 7T < Ce~ T/2 coshs, |cj t | 7t < Cr 3 e~ 3T/2 . 

The modifications needed for the asymptotically Euclidean and asymptotically hyper- 
boloidal cases are discussed in §1?] and statements of this theorem in these settings are 
given there, in Theorems |] and ^ 

Notice that we have not stated the precise regularity of the solutions, but as indicated 
earlier, we work in various classes of Holder spaces; in particular, we also obtain existence 
of initial data sets (£,7, II) with only finite regularity. 

Some estimates of the geometry of the solutions are provided in §8; these should be 
of interest when considering the spacetime evolutions of these initial data sets. Finally, 
§9 contains an informal discussion of various concrete examples of initial data sets, many 
of which are new, which our construction gives. 

2 The approximate solution 

In this section we construct the family of approximate solutions ( 7 t, [Mr, V't)- Let rj 
denote the geodesic distance to pj (relative to 7) in the ball Bj. By Gauss' lemma we 
have 

tIbj = dr 2 + r 2 h j (r j ), 

where hj is a family of metrics on S 2 , smooth down to rj = and such that hj(0) is 
the standard round metric on the sphere. The standard way to obtain a metric with 
asymptotically cylindrical ends is to multiply 7 by rj 2 in each of these balls. Recalling 
the convention of (|j) that the conformal factor is written as ip 4 ", let us fix a function ip c 
which is strictly positive on £*, which equals one on anc ^ which equals r 1 - 2 on By 
Now define 

7c = C 4 7- (8) 

This metric agrees with 7 on and has asymptotically cylindrical ends in a neighbor- 
hood of each pj . 
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Figure 1: Conformally blowing up (£,7) to form (S*,7 C ) with two asymptotically 

cylindrical ends. 

In fact, setting tj = — log rj, we see that 

lc\ B] =dt) + h j {e- t i). 

Recalling that hj is smooth in rj, we can rewrite this as 

dtj + d0 2 + e~^kj, (9) 

where d9 2 is the standard round metric on S 2 and the term kj in the remainder is 
a symmetric two-tensor which is smooth in rj = e~ tj . Note that if the metric 7 is 
conformally flat in Bj then the definition of ip c may be easily modified to make the 
metric j c exactly cylindrical in Bj. Note also, for later reference, that 

-logi?< tj < 00. (10) 

Next, following @, define 

Mc = V>cM; (11) 
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this is defined on S* and is transverse-traceless with respect to j c . Relative to the initial 
data set (7 c ,/i c ) on £*, if; c satisfies the Lichnerowicz equation 

A c ^ c -\R c i> c + \ |/i c | 2 " ^ r 2 ^ c = 0. (12) 

In this equation, the expressions A c , R c and | • | 2 are all computed in terms of the metric 
7 C . Thus (7c, /j, c ,ifj c ) is an initial data set on £*. 

Let A = — log i? and let T be a large parameter. Truncate the manifold S* by omitting 
the regions where tj > A+T. We may form a new smooth manifold by identifying the two 
finite cylindrical tubes {(tj,6) : A < tj < A + T} via the map (ti,9) —> (T — 1\, -9). Let 
us call the resulting manifold St and let the now-identified cylindrical tube be denoted 
Ct- (Of course, each £y is diffeomorphic to S, the manifold obtained by adding a handle 
to S, but it is convenient to keep track of the dependence of T explicitly.) We introduce 
a new linear coordinate s on Ct by setting 

s = ti - A - T/2 = -t 2 + A + T/2. 

Thus Ct is parametrized by the coordinates (s,9) <E [—T/2, T/2] x S 2 . 




Figure 2: The manifold (St,7t)- 

To continue we must define the family of metrics fr, 'almost-transverse-traceless' 
tensors jjlt on Sy, and conformal factors ipT, which we do in turn. 

The definition of 7t is most transparent when 7 is conformally flat in the balls Bj . In 
this case, using the modifications of ip c alluded to above, in the cylindrical coordinates 
(tj,9), we have 7 C = dt 2 + d9 2 , and thus the map identifying the two cylindrical pieces 
to Ct is an isometry and 77 1 is thereby well-defined. In the general case, we cover T,t 
by two open sets U\ and U2, the intersection of which consists of two components, one 
disjoint from B\ U B2 and the other equal to {(s, 9) : — 1 < s < 1}. Choose a partition of 
unity {xi>X2} subordinate to this cover and let jj and Hj denote the restrictions of 7 C 
and n c to Uj. Then, with the obvious abuse of notation, define 

IT = Xl7l + X272, MT = XlMl + X2^2- 

Notice that jt = 1c and \it = l~i c away from the middle Q = [— 1 , 1] x S 2 of the neck 
region, where ti,t 2 ~ T/2. 
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Figure 3: Defining the metric 7t = Xi7i + X272, and the approximate TT-tensor 

= XlA*l +X2/"2- 



Although it is not obvious why at this stage, we need to define ipx somewhat dif- 
ferently. This time choose a covering {Z^i,^} of S such that Ui(~)U2 C S^ and with 
Pj £ ^j. Also choose nonnegative functions %j £ C^iUj) with Xi + X2 = 1 ° n SJj, and 
such that when restricted to Bj, Xj = 1 for tj < A + T — 1 and Xj = for tj > .A + T. 
Writing the restriction of ipc to t/j as tpj , we set 

V>T = XlV'l + X2^2- 

Then, based on the construction of S^ from the conformal blow up of (S, 7), we define 
V't on in the obvious way. Note that ipT is exactly equal to 1 away from the cylinder, 
and is equal to tpi + ^2 throughout most of the cylinder except at the two ends. At these 
ends, near the junctions of the cylinder with S R , one of the two functions tpi or ^2 is 
roughly of unit size and the other is exponentially small; in defining ifr, the exponentially 
small summand is cut off to be zero near these junctures. 




Figure 4: The approximate solution to the Lichnerowicz equation, ipT = XiV'i + X2^2- 

In summary, we have defined 7^ and by using cutoffs supported on the middle of 
the cylinder Ct, but have cut off either of the tpj at the ends of the cylinder. The reasons 
for doing this are that j c decays to the product metric like e~ tj and fi c is extremely small 
in the middle of Ct, so the error terms from these are of order e _T//2 . On the other 
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hand, ip c only decays like e~ tj ^ 2 along each end, and so if we were to cut this function 
off in the middle, it would produce an error term with size approximately e _T//4 , which 
is unacceptably large. We carry out these estimates in more detail in §|j. 



3 Transverse-traceless tensors on 

We now undertake the first of several steps intended to perturb the approximate solution 
(TTjMTjV't) constructed in the last section to an initial data set, i.e. an exact solution, 
when T is sufficiently large. Our goal in this section is to modify the tensor so that 
it becomes transverse-traceless with respect to 7t- 

Before starting, note that it is clear from the definitions that (i? is transverse-traceless 
away from the center of the neck region Q in the middle of the cylinder Ct- Furthermore, 
by choosing the cutoff functions carefully, we may also assume that \xt is trace-free on 
all of Tit] however, its divergence is almost surely nonzero, but at least is supported in 
the annulus Q. 

There is a well-known procedure for producing a correction term to kill this diver- 
gence. We recall this now, and refer to |l5| or [33] for more details. Let 7 be a metric 



and X a vector field on the three-manifold S, and write 

VX= l -C xl -^{d\vX) 1 . 

Equivalently, in local coordinates, 

(VX) jk = - {Xj, k + X k .j) - ^div (X) ljk . 

The first order operator T> = 2? 7 defined by this expression maps vector fields to trace- 
free symmetric (0, 2) tensors and has the property that T>X = if and only if X is a 
conformal Killing field. The formal adjoint of T> is T>* = — div , and the second order 
operator L = T>* o T> is self-adjoint, nonnegative and elliptic. This operator is often 
referred to as the 'vector Laplacian'. 

Now let 7 = 7T) with T> and L the associated operators. Let W be the vector field 
associated to div \it using the natural duality between vector fields and 1-forms. Suppose 
that we can solve the equation LX = W. Then writing ot = L)X we compute that 

div Or — o"t) = W — divVX = W — LX = 0, 

and so fir — ot is divergence- free. Since \it and or are both trace-free, we have produced 
the desired correction term. 

To implement this strategy properly, we must show that when T is large enough, 
the operator L is invertible, and moreover that the solution X, and hence T>X, is much 
smaller than \it, so that ot can honestly be regarded as a small perturbation. Note that 
since W is a divergence and hence orthogonal to the cokernel of L, we can always solve 
LX = W (but without an estimate on the solution). The main issue is to prove that the 
inverse G to L on the manifold Sy, when acting between appropriate function spaces, 
has norm bounded at worst by some polynomial in T. (This is sufficient because the 
error term div (i? decays exponentially.) 
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When 7 is conformally flat in the balls Bj, it is possible to choose the conformal 
factor tp used to define j c slightly differently so that is the product metric on the 
cylinder C T = [-T/2,T/2] xS 2 . We proceed in this case by an analysis of the operators 
T> and L on the complete cylinder R x S 2 as well as on the finite cylindrical piece Ct , 
using separation of variables and explicit calculations with the resulting ODEs, to reach 
the conclusion that the inverse G to L exists and has norm (as an operator on L 2 (Et)) 
blowing up no faster than T 2 . When 7 is not conformally flat in these balls, an additional 
perturbation argument is required to reach the same conclusion. The basic goal in these 
arguments is to show that the lowest Dirichlet eigenvalue of L on Ct is of order T~ 2 . 
Afterwards, the same conclusion is reached for the operator L on Sy. 



3.1 The conformal Killing operator on the cylinder 
3.1.1 The decomposition 

We first derive the explicit form of L on the cylinder Kx5 2 with product metric ds 2 + h, 
where s is the linear coordinate on R and h is the standard round metric on S 2 . Any 
vector field X on the cylinder may be written 

X = fd s + Y(s) 

where / is a function on the cylinder and Y(s) is a vector field on the sphere {s} x S 2 . 
Write VX = S(X) — ^divX7, and choose any local coordinates 6 = (61,62) on § 2 . Then 
for i,j = 1,2, 

S(f d s ) 00 = d s f, S(fd s ) 0i = S(fd s )ij = 

and div(/<9 s ) = d s f. Therefore, exhibited as a matrix, 

jd s f \V e f 



V(fd s ) 



\V e f -h(d s f)h 



where Vg represents the gradient of / on the sphere. Hence, with the convention Aq 
-div Vo, 



L(fd s ) = -div(D(/aO) 
Next, 



-ld 2 f+\A e f \ /-|a 2 / + lA e/ 
-kds(V e f) + lV e (d s f) J \ -$d a (Vef) 



S(Y) 00 = 0, S(Y) 0i = -d s (Yi), S(Y) ZJ = SeWij, 



1 

where by definition the conformal Killing operator T>h on the sphere decomposes as 
V h Y = S e (Y) - \{d\-v e Y)h. Using also that div(Y) = div e (Y), we get 



V{Y) = 



-idiv fl (y) \d s Y 

\d s Y V (Y) + ldiv (Y)h 
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and then that 

/ i9 s (div e (y)) - idiv e (d s y) 

L(Y) = -div(D(F)) = 

V -\d 2 s Y - &v e {V e (Y)) - iV e div e (y) 
-\d s {div e {Y)) 

-\8 2 S Y + L e {Y)-\V e d\v e {Y) 

It is possible to express this operator in a somewhat more familiar form. To do this, 
we identify vector fields and 1-forms on S 2 in the usual way using the metric h. Thus, 
Vg and dive correspond to dg and —Sg, respectively. We still let Y denote the 1-form 
associated to the vector field Y, hopefully without causing undue confusion. Using these 
identifications, 



L(fds + Y(s)) 



-ld 2 f + ^gf + \d s {5gY) 

-\8 2 Y + Lg(Y) + \dg5g{Y) - \d s (dgf) 
This expression may be simplified using two separate Weitzenbock formulae, 

Lg = ^V*V-^ : and Ag = V*V + 1, 

both of which use the fact that the curvature tensor on S 2 is constant. The first of these 
relates Lg to the 'rough Laplacian' V^Vg (cf. [33]) and the second is the more familiar 



one relating this rough Laplacian to the Hodge Laplacian Ag = dgSg + 5gdg. These 
combine to give 

Lg = ^Ag - 1, 

and thus we finally arrive at the expression 





\d 2 + ^dg 5g + ^Sg . 




(13) 



3.1.2 Separation of variables 



To analyze this operator L on the cylinder in more detail, we introduce the eigenfunction 
expansion on the S 2 factor. Let be an orthonormal set of eigenfunctions for the 

scalar Laplacian on S 2 , so that Agcf)^ = Xk4>k an d each A/% = j(j + 1) for some j £ N. 
If * is the Hodge star on S 2 , then {X k l ^ 2 dg4>k-, \ k l ^ 2 5g * 4>k} is an orthonormal basis of 
eigenfunctions for the Hodge Laplacian on 1-forms. Set (j) = <j)k and A = Xk with A& > 0, 
for some k, and consider a 1-form of the form 



X = f(s, 8)ds + Y(s) = u{s)<j)ds + v(s)^ + w(s)-^jt, 

v A v A 
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or equivalently, the associated dual vector field. Then 

LX = {-\u" + ±u+ ^ d>ds + (-\v» + (\\ - l)v - ^u'j 

+ {-2 W +i 2 X - 1)W )^T> 
and so L acts on the column vector (u, v, w) 1 by 

-| o \ / ^o\ / I o 

-\ 9 S 2 + o 9.+ |A-1 

-\ j \ o 0/ V ±A-1 

There is one extra case, when A = Ao = 0; then X = u(s) ds and LX = —^u"(s)ds, so 
we write 



L = -Idl 



3 

When j > 1, L,- uncouples as Lj = L'- © L" where 



\ n o . f ^ 



■i-V o 4T sT l-f o r s "Vo 1A-1 



acts on (u, u)* and 



acts on w. Thus altogether, 



^' = -^ + (^-1) 



L = L © 0(L^ © L'j). 



3.1.3 Temperate solutions 

The first application of these formulas concerns the solutions of LX = on the cylinder 
which do not grow exponentially in either direction. 

Proposition 1 The space of temperate solutions to LX = on the cylinder R x S 2 is 
8 -dimensional and is spanned by the 2- dimensional family of 6 -independent radial vector 
fields (a + bs)d s , a,b Gt, and the 6- dimensional family of vector fields which are dual to 
the 1-forms (q + dis)5g * <pi, Ci, di G R, where 4>i = xi\ S 2, i = 1, 2, 3. 

Proof: It is straightforward to check that these bounded or linearly growing vector 
fields are all in the nullspace of L. In fact, using the eigendecomposition of L above, 
Lqu = = implies u(s) = a + bs, and this gives the first family. Next, for j = 1 

(i.e. A = 2), L"w = —\w" = implies w(s) = c + ds; the eigenfunctions here are the 
restrictions of linear functions on M 3 to S 2 , and so this gives the other family. It remains 
to show that these are the only temperate solutions. 
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It is clear that when j > 2, hence A > 6, all solutions of L"-w = grow exponentially 
in one direction or the other. On the other hand, consider the homogeneous system V- 
when j > 1, which is 

2 „ A v A , 
_„»__„__„< = 

+ = 0, 

where A = j(j + 1). This is a constant coefficient system, and solutions have the form 
" 1 e ps z 1 , or se ps z x + e ps z 2 , z 1 ,z 2 £C 2 , 



where p is a root of the indicial equation 



2^2 _ A -^Ap 

del [ " | =i(p 4 + 2(l-A)p 2 + (A 2 -|A))-n. 

„ 1 J2 (2- 



^-(|A-1) 



3 



Solutions are temperate only when Re p = 0, but it is easy to check that this polynomial 
for p has no root with purely imaginary part, since A = j(j + 1), j > 1. The proof is 
complete. □ 

Remark 2 T/ie space of conformal Killing fields on S 3 is 10- dimensional, and is gener- 
ated by the 6- dimensional family of rotations and the 4-dimensional family of dilations. 
Amongst all of these, only the vector fields vanishing at the two antipodal points ( cor- 
responding to the infinite ends on the cylinder) give rise to temperate conformal Killing 
fields on R x S 2 . Of course, these are just the generators for the 3-dimensional family 
of rotations of the orthogonal S 2 (i.e. corresponding to 5q * (pi, i = 1,2,3, as above) and 
the dilation fixing these antipodal points (d s ). 



3.1.4 The lowest eigenvalue of L on [-T/2,T/2] x S 2 

The next application concerns the lowest eigenvalue of the operator L with Dirichlet 
boundary conditions X = on the finite piece of the cylinder Ct = [— T/2,T/2] x S 2 . 

Proposition 2 Let Xq = \o(T) denote the first Dirichlet eigenvalue for L on Ct = 
[-T/2,T/2] x S 2 . Then 

for some constant C independent of T . 

Proof: This estimate is clearly sharp for the scalar ordinary differential operators Lq = 
— 1<9 2 and L" = —\d 2 . On the other hand, when j > 2, the lowest eigenvalue of L" 
converges to + 1) — 1 as T —> oo, hence is bounded away from zero. 



14 



We finally show that the lowest eigenvalue for L'- is also bounded away from zero 
when j > 1. In fact, 

i; (")•(")) = C ^ ,)2 + ^ + ^ iuv ' - u ' v) + ^ + & - i>2 * 



T/2 2 ( , V\ Y if, V\ Y 35 x 2 ,63 x 1N 2 , 
_ T/2 3^ -- W J + 2^ + T U J + ~2 XU +( 96 A " 1)W ^ 

Since A > 2, this is bounded below by \ Jq(u 2 + v 2 ) ds, and the proposition is proved. □ 
3.1.5 Growth estimates on Ct 

To conclude this analysis of L on the cylinder, we show that if X is an eigenfunction 
with 'very small' eigenvalue, then its size in the interior of Ct is controlled by its size on 
the ends of Ct- The next result does this directly for the product metric. After that we 
prove a more technical estimate for the inhomogeneous eigenvalue problem which will be 
used later to handle the general case by perturbation. 

Proposition 3 Suppose that LX = fiX, where L is computed with respect to the product 
metric and \i < |Ao(T). Suppose also that 

[ (\x(-t/2, e)\ 2 + \x(t/2, e)\ 2 ) de < c x . 

Js 2 

Then for any a G [-T/2 + 1, T/2 - I], 

ra+l p 

/ / \X{s,6)\ 2 dsdO < C 2 , 

Ja-l JS 2 

where C2 depends only on C\, but not on T or a, and where C2 — ► as C\ — > 0. 

Proof: It suffices to check this estimate on each of the components in the decomposition 
of L. For the scalar operators Lq and L" this estimate may be verified using the explicit 
solutions of the equation (though this can easily be done by other standard and less 
computational methods too). So as before, it suffices to focus on the system Lj, j > 1. 
Write 



L' 3 = -Ad 2 + Bd s + C 



where 



a _ f 2/3 \ / VA/6 \ _ / A/2 

V 1/2 ; ' V -VA/6 ; ' V ((2/3)A - 1) 

and, as usual, A = j(j + 1). 

Now write X = (u,v) 1 and define f(s) = (AX,X). Then 

l -d 2 s f(s) = {AX",X) + {AX',X'). 
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Using the equation L'-X = fiX and the skew-symmetry of B, this becomes 

\d 2 J(s) = -(X',BX) + (CX,X) + (AX',X') - n(X,X), 
which can then be rewritten as 

l -d 2 J{s) = \AV*X' - \a^I 2 BX\ 2 + (DX,X), 

where 

But D is positive definite when A > 2 and T is large enough that subtracting fi does not 
destroy positivity, and so we conclude that f"(s) > 0, and hence that the L 2 norm of X 
at s = ±T/2 controls the L 2 norm of X over any strip a — 1 < s < a + 1. □ 

Our second result here is more specialized, and is the precise additional ingredient 
required later to handle the nonproduct case. We let || • ||o,a,<r denote the C 0,a norm over 
the set [a, a + 1] x5 2 . 

Proposition 4 Let L be computed relative to the product metric on Ct, and suppose 
LX = fiX + F, where X(±T/2) = and fiT 2 < uq for some sufficiently small constant 
uq independent of T . Suppose also that for all \s\ < T/2, 

\\F(s,6)\\ 0:a:a <f{a)e' T l 2 cosh a, 

where |/(c)| < 1 for all a. Recalling that F implicitly depends on T, we also assume that 
for any e > and A > 0, there exists a Tq such that when T > Tq we have |/(<r)| < e 
when T/2 — A < \a\ < T/2. Then, for any r) > we have \X\ < r] on all of Ct, when T 
is large enough. 

In other words, if X solves this inhomogeneous eigenvalue problem on Ct, with 
vanishing Dirichlet data, if /jl is much smaller than 1/T 2 , and if the forcing term F is 
uniformly small near the ends of Ct and bounded by e~ T l 2 cosh s on the middle portion 
of this cylinder, uniformly as T — > oo, then X is uniformly small on all of Ct- 

Proof: We separate variables and write the equation as LjXj = Fj for j > 0. For 
j > this equation splits further into LjX'j = Fj and L'-X'- = F" , as in the previous 
subsection. None of the primes here correspond to derivatives. 

We have shown that the operators L'j with j > 1 or L" with j > 2 have spectrum 
bounded below by a positive constant as T — > oo. For any of the corresponding equations, 
if the conclusion were false then we would be able to produce, by the usual arguments, 
a nontrivial bounded function X' or X" such that L'-X' = or L"X" = 0, respectively, 
defined on the complete cylinder M x S 2 . This contradicts the strictly positive lower 
bound for the spectrum. 

For the remaining components, it suffices to prove the assertion for the scalar problem 
d 2 u + a 2 u = F, where aT is sufficiently small. For notational convenience we shall use 
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the shifted variable t = s + T/2 (so that F(t) = f(t)e~ T/2 cosh(t - T/2), where this 
function / is a translate of the one above. We may write the solution explicitly: 



sin a(T — t) f l sin ar„, , , sin at /" T sin a(T — r) „. . , 
u(t = / T 1 / ~^F(t) dr + — - / . V T ' F(r) dr. 

Because aT is so small, we may approximate each of the sines with their first order 
Taylor approximations; using the estimate for F too we obtain 



\u{t)\ < C 



J\ cosh(r - T/2)e- T / 2 /(r) dr 
+t £ (l - cosh(r - T/2)e~ T l 2 f{r) dr^ . 



Finally, we must use the separate upper bounds for f(r) in the various regions < r < A, 
A < t < T — A and T — A < r < T. The main fact used in this estimation is that the 
function h(t) = t(l — t/T)\ sinh(i — T/2)|e _T//2 is uniformly bounded independently of T, 
attains its maximum near t = 1 and t = T — 1, and for any 5 > there exists an A > 
such that h(t) < 5, uniformly in T, when A < t < T — A. Together with the bounds for 
/, we can derive that u is uniformly small. Details are lengthy but straightforward, and 
are left to the reader. □ 



3.2 The lowest eigenvalue of L on E T 

Theorem 2 The lowest eigenvalue Ao = Ao(T) for L on satisfies Xq > CT~ 2 for 
some constant C > which is independent ofT. 

Proof of Theorem ^| when the restriction of jt to Ct is the product metric: 

If this result were false, it would be possible to find a sequence of values Tj — ► oo such 
that for the operators Lt 3 associated to the metrics 7^., we would have Xo(Tj)T 2 — ► 0. 
We now show this leads to a contradiction. 

Denote the sequence of metrics and operators by jj and Lj, respectively. Let Xj 
be any one of the eigenfunctions for Lj corresponding to the eigenspace with eigenvalue 
Xo(Tj). By rescaling, assume that sup ST \Xj\ = 1. 

The coefficients of Lj converge uniformly on any compact subset of £* to the coeffi- 
cients of the operator L c corresponding to the metric j c . Since \Xj \ is uniformly bounded, 
we may extract a subsequence Xj> which converges uniformly on compact sets in S* to 
a limit X. This vector field is a solution of L C X = and satisfies \X\ < 1. 

Our first claim is that X is nontrivial. To see this, relabel the subsequence as Xj 
again. Suppose that the supremum of \Xj\ is attained at some point qj € Ey.. If qj lies in 
E^ \ Ct 3 for infinitely many j, then again possibly passing to a subsequence qj — > q G E* . 
Since = 1, we get = 1, and so X ^= 0. On the other hand, suppose that 

qj G Ct for all but finitely many j. Note that by Schauder theory, \Xj \ remains bounded 
away from zero on a strip [oj — r], dj + 77] x S 2 C CV of width not depending on about 
qj. Furthermore, Xq(Tj) < ^T 2 for j large enough, and hence Proposition ||| (and elliptic 
theory) implies that |Xj(g'-)| > a > for some point q'- £ dC^ and some a not depending 
on j. Therefore in this case too we may conclude that the limit X on E*. 
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By the analysis of §3.1, if X is any solution of L C X = on £* (with respect to the 
metric j c ) which is bounded along the cylindrical end, then it can be written as a sum 
X = Xqo + Z, where extends to a bounded solution of LX^ = on the entire 
cylinder R x S 2 (with product metric) and Z decays exponentially as t —* oo. 

Our second claim is that, assuming the nondegeneracy of S so that there are no 
conformal Killing fields vanishing at both points pj, then any bounded vector field X on 
£* which satisfies L C X = necessarily vanishes identically. To see this, observe from our 
earlier result about elements in the nullspace of L c which are bounded on the cylinder, 
all such solutions are also annihilated by T>. (This is not true if we also admit the linearly 
growing solutions.) Thus the term X^ in the decomposition of X above is annihilated 
by T>, and hence T>X = DZ decays exponentially along the cylindrical ends. Therefore, 
no boundary contribution appears in the integration by parts 

0= (L C X,X) = (D*VX,X) = \ \VX\\ 2 , 

so that T>X = 0. Thus X is a conformal Killing field on £* which is bounded with 
respect to the metric 7 C ; changing back to the old polar variables rj = e~ tj and 6 on 
each ball Bj shows that X is bounded with respect to the original metric 7 on X, and in 
fact vanishes at both points pj. It is therefore in the (distributional) nullspace of L on 
S, and hence by elliptic regularity extends smoothly to a solution on all of S. It would 
thus be a nontrivial conformal Killing field on the manifold S, vanishing at both p\ and 
P2, and we have assumed that such a field does not exist. This completes the proof. □ 

Proof of Theorem [| in general: We may follow the outline of the preceding proof 
at almost every point. Thus we suppose the result false and choose a sequence of eigen- 
functions Xj with corresponding very small eigenvalues fij, for some sequence Tj — > 00. 
Supposing that \Xj\ < 1, we again attempt to pass to a limit and obtain a contradiction. 
If \Xj\ remains bounded away from zero on some fixed compact set near the ends, then 
we argue as before. The difficulty, however, is when \Xj\ remains bounded away from 
zero in some portion of the cylinder Ct^ which recedes from either end. To handle this, 
we would like to use the same sort of argument as in the proof for the product case. 
Namely, we would like to assert that if \Xj\ tends to zero near the ends of Ct^ then it 
must tend to zero uniformly along all of this cylindrical segment. We shall now prove 
this using Proposition ||. 

The metric 77 1 differs from the product metric on Ct by a term of size e _T / 2 cosh s, and 
therefore the operator L computed relative to 77 1 differs from the corresponding operator 
Lq, computed relative to the product metric, by an error term Et with coefficients of 
this same size. We may rewrite the equation for Xj as 

L Xj = fijXj + Fj, 

where Fj = —E^Xj. Now suppose that we are in the difficult situation: that \Xj\ is 
uniformly small on any fixed neighborhood of the ends of the cylinder. Then it is clear 
that Fj satisfies precisely the estimates in the statement of Proposition |j. The solution 
Xj may be decomposed into a sum Uj + V}, where Uj satisfies the homogeneous problem 
Lollj = HjUj and Uj(±T/2) = Xj(±T/2), while Vj satisfies the inhomogeneous problem 
with vanishing boundary conditions. Clearly both of these summands are uniformly 
small near the ends of Cx r The uniform smallness of |t7j| on the entire cylinder now 
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follows from Proposition ||, while the uniform smallness of | Vj | follows from Proposition ^. 
Therefore it is impossible that Xj is large in the middle of Ct s but small at the ends. 
This completes the proof. □ 

3.3 Estimates on Holder spaces 

We have shown that the lowest eigenvalue of the operator L on is bounded below 
by CT~ 2 , uniformly as T — > oo. Hence the norm of G = L , as a bounded operator 
on L 2 (E T ), blows up no faster than C'T 2 . However, it is more convenient for us to use 
Holder spaces because of the nonlinear nature of the problem, and we show now that 
this L 2 estimate may be converted to a slightly weaker estimate on appropriately defined 
Holder spaces. 

Definition 2 Consider the cover of given by the sets S^/ 2 and Ct, and for any 
vector field X on St, let X' and X" denote the restrictions of X to these two subsets, 
respectively. Also, for -T/2 + 1 < a < T/2 - 1, define N a = [a - 1, a + 1] x S 2 C C T - 
Then for any k £ N and a £ (0, 1), define 

\\X\\k, a = \\X'\\k, a + SUp ||X ||fe, a ,JV a - 

-T/2+l<a<T/2-l 

The local Holder norms here are computed with respect to the metric 7t- 

Corollary 1 Let W £ C k ' a {T, T ) and suppose that X is the unique solution to LX = W , 
as provided by Theorem ^. Then for some constant C independent of W and T, 

\\X\\k+2,a < CT 3 ||VF|| fcja . 



Proof: First, by Theorem |, \\X\\ L 2 < CT 2 \\W\\ L 2. Next, since the volume of S T is 
bounded by CT, we have ||W||x,2 < T||W||o,q. In addition, local Schauder estimates for 
L c give ||-X"||2,a < C(||W||o,a + II-^IIl 2 )- Putting these together yields 

||^|| 2 ,a<cr 3 ||H/|| 0ia , 

which is the desired estimate when k = 0. The estimate when k > requires another 
application of the local Schauder estimates. □ 

3.4 Correcting fj,x to be transverse-traceless 

Recall now that fix is constructed by patching together fi c using cutoff functions in the 
center Q of Ct, and the metric 7t is constructed from 7 C the same way. The definition 
(|l~l| ) implies that on Q 

Imt| 7t < cv4m 7 ( 14 ) 

with equality (for C = 1) on \ Q. Since \i c is divergence-free with respect to j c , the 
'error term' div 7T //-r arises from two sources: the cutoff functions used to define /it and 
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those used to define jx- This error term is clearly supported in the center Q of Cx, and 
it is not hard to see that 

||div 7T /i T || fc _ lja < C'\\fj, T \\k,a < C"\\ip T \\ k>a , 



the last inequality following from (14) since |//L < C. Recalling that ifo ~ e~ 1 ^ 2 and 
each tj ~ T/2 on Q, and letting W be the vector field dual to div 7T /iT, then we conclude 
that 

\\W\\ k+1 , a <Ce- 3T / 2 . (15) 

Proposition 5 There exists a tensor ax on Y<x with 

\\o- T \\k,a<CT 3 e- 3T / 2 (16) 

and such that fix = fJ>T — o~t is transverse-traceless with respect to jx- 

Proof: As indicated in the beginning of this section, we first solve the equation LX = W 
on Sy, with W = div 7T ^, and then set ax = T>X. The operators T> and L here are 
both computed with respect to jx- The estimates for gt and fix follow from ([To]) . □ 

4 Estimates on the approximate solution 

Our overall goal in this paper is to modify the approximate initial data set (7^, (j,t,iPt), 
defined on the manifold Et, to a genuine initial data set (77 1 , At> ipr), so that fix is 
transverse-traceless with respect to 77 1 and i()t satisfies the Lichnerowicz equation relative 
to 77 1 and fix- Thus far we have accomplished the first of these tasks, and so what remains 
is to show how to modify the conformal factor ipx to ipx = V'T + f]T so that ipx satisfies 
the Lichnerowicz equation. There are three separate steps involved in this. In the first 
we must calculate the 'error term' Ex, which measures the deviation of ipx from solving 
the Lichnerowicz equation. After that we must analyze the mapping properties of the 
linearization of this equation, calculated at ipx, uniformly as T — > 00. Finally, these 
are put together to carry out the contraction mapping argument to produce the exact 
solution ipx- These steps are taken up in this and the next two sections, successively. 
To begin then, denote by Mx the Lichnerowicz operator with respect to {^Xii^r)'- 

N T {i>) = A T ^ - \Rx^ + \\M 2 ^~ 7 - ^V- (17) 

The Laplacian and scalar curvature and norm squared of fix are all computed with respect 
to jx here. We define the error term 

E T =Mx(^x), (18) 
which measures the deviation of ipx from being an exact solution. 

Proposition 6 For every k and a there exists a constant C which does not depend on 
T such that ||£r||fca < Ce~ T / 2 . 
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The proof of this proposition involves a careful accounting of the various terms in 

We begin by noting that away from the neck region Ct, the estimate is nearly trivial. 
In fact, on T,* R , jt = 7o ipT = V'c and P-T = He + 0(T 3 e~ 3T l 2 ). Inserting these into the 
Lichnerowicz equation, and using that (-j c , /j> c , ip c ) is an exact solution, we have that 

\E T \ < CT 3 e- 3T / 2 < Ce~ T / 2 , 

along with all its derivatives, as claimed. 

We divide the tube Ct up into three regions: the center Q = [—1, 1] x S 2 , and the 
regions to the left and right of this, [-T/2,-1] x S 2 = C T l) and [l,T/2] x S 2 = C T 2) . 
We also recall the notation used in the construction of the approximate solution that jj, 
\ij and tjjj denote the restrictions of 7 C , \x c and ip c to the balls Bj, and hence, by a small 
abuse of notation, to the tube Ct as well. Thus (71, (j,%,ipi) represents data 'coming in 

from the left', while (72, A*2, ^2) represents data 'coming in from the right'. Recall also 

(i) 

that 7t and [it agree with jj and fj,j on Cj, , whereas i/jt is a sum of two terms, one 
much larger than the other, in each of these side regions. 

We begin by listing some trivial estimates which hold in all of Ct, and which follow 
from the construction of the approximate solution and from (16): 

a) 7T = ds 2 + h + 0(e~ T/2 coshs), 

b) ~ 2e- T / 4 cosh(s/2), 

c) |^r| ~ ip T < Ce- 3T / 2 (cosh(s/2)) 6 , 

d) |At-Mt| <CT 3 e~ 3T / 2 . 

We shall also repeatedly use the fact that 



A 7j ^ - 7^(7;)^ + M% 7 ~ = 0- (19) 



We begin by estimating Et in Q. The main observation concerning this estimate is 
that we have a decomposition 

i>T = i>\ + 1P2 

valid in Q. In addition, it follows from property (a) above and from the definition of 77 1 
that 

A T -Ii? r = A 7j -^ 7j +0(e- T / 2 ) 

for j = 1 or 2. On Q, we may therefore expand (A^ — \Rt)^>t as a sum of two terms. 
Since on Q, ipj = 0(e~ T ^), we see from ( |l9| ) that we have 

(At - \Rt)^t = X>~N V + ^V?) + 0(e" 3T/4 )- 
From this, using the estimate that \^j\ 2 = 0{e~ 3T ) in Q, we easily conclude that 

(A T - \r t )4>t = 0{eT 3T l A ). 
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As for the other terms in the Lichnerowicz equation, using properties (b), (c) and (d) 
above we easily get suitable estimates for them in this region as well. In particular, we 
have \fi T \ 2 ^T 7 = 0(T 6 e' 5T / A ) and t 2 ^% = 0{e- 5T / 4 ). Hence in Q, we obtain \E T \ = 
0(e~ 3T / 4 ), which is even better than the estimate we have stated. 

The estimates in C^p and are nearly identical to one another, so it suffices 

(2) 

to concentrate on just one of these, let us say in C T to be concrete. In this region, 
4>T = 4>2 + xi^ii where x is a function cutting tpi off to zero at the far end of Ct, where 
s « T/2. In this region 

^ = e - T /V s / 2 + 0(e- T / 4 e" 3s / 2 ), ^ ~ e^'V' 2 , 

and in addition 

A T - \r t = A - -Rq + 0{e~ T / 2 cosh s). 

o o 

Therefore 

(A T - Ir t )Mi) = e- 3T/4 cosh S e-*/ 2 + 0{e~ T l 2 ) = 0{e~ T ' 2 ). 
8 

Next, ipT = ^2(1 + X&i/fa) and hence 

^ = ^ + 0(e 7T ' A e~ 9s / 2 ), 

and 

V- T = ^ 5 + 0(e- 5T / 4 e 3s / 2 ). 

Furthermore, we have here 

\t T \ 2 = \V2\ 2 + 0(T*e- 3T / 2 ) 

Putting these estimates all together yields \Et\ < Ce~ T / 2 on cf^\ With a similar 

estimate on C^ 1 "* we conclude that \Et\ < Ce _T / 2 on all of Ct- The estimates on the 
derivatives of Et follow similarly. 

5 The linearization of the Lichnerowicz operator 

Our main goal now is to find a correction term tjt for the approximate solution tpx so that 
the sum ip? + Vt solves the Lichnerowicz equation Mt^t + Vt) = on St- This will be 
done in the next section using a contraction mapping argument. In preparation for this 
we establish in this section good estimates for the mapping properties of the linearization 
Ct of this operator Mt about the approximate solution ipT] as usual, keeping track of 
behavior as T — ► 00. 

We first record that the linearization C of the Lichnerowicz operator on £ with respect 
to the initial data set (7, fx) is given by 

£ = A 7 -^(i?( 7 ) + 7H 2 + yT 2 ), 
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but using the fact that the pair (7, /i) satisfies the constraint equations, we can rewrite 
this as 

C = A 7 - (\fj,\* + lr 2 ). (20) 



7 V ^ 3 
By standard elliptic theory (for E compact), 

C : C fc+2 ' a (S) -> C fc ' Q (S) 

is a Predholm mapping of index zero for every k £ N and a £ (0, 1). Furthermore, since 
we are assuming that II ^ 0, which is equivalent to r 7^ or fi ^ 0, the term of order zero 
here is nonpositive and not identically zero and so the maximum principle implies that 
there are no nontrivial solutions of C<j) = 0. Consequently the cokernel is also trivial, and 



we see that (20) is always an isomorphism. 

Our first goal in this section is to establish that the same property is true for the 
linearization Ct on Ey, which is given by 



C T = A 7T - i ( R{ 1T ) + 7 |/i T | 2 ^ 8 + y r 1 < i ). (2.1 

This operator is elliptic and Fredholm of index zero acting on the Holder spaces intro- 
duced in §3.3; unfortunately, the term of order zero here is no longer necessarily nonposi- 
tive since it can no longer be rewritten using the constraint equations, and we do not have 
control on the sign of R{^t) near the boundary of T,* R . Therefore, the triviality of the 
nullspace of Ct when T is large enough requires further argument. We must also show 
that the norm of the inverse operator is uniformly bounded as T — > 00. This last issue is 
complicated by the following considerations. We shall be using this solution operator in 
an iterative scheme to produce the correction term t\t for the approximate solution ipx- 
Because of the term t/j~ 7 in the Lichnerowicz equation, it is quite important that r\T be 
substantially smaller than ipT- Furthermore, we also desire that r]T vanish exponentially 
as T — > 00 on the 'body' ~S* R , away from the scene of the surgery. To ensure that both of 
these requirements are met in the iteration scheme, we establish precise estimates for the 
solution operator for Ct- This necessitates using slight alterations of the Holder spaces 
which incorporate a weight factor along the neck. Thus we first define these weighted 
Holder spaces and then take up the matters of the existence and uniformity of the inverse 
of Ct on these spaces. 

Definition 3 Let wt be an everywhere positive smooth function on Y<t which equals 
e —T/4 CO sh(s/2) on Ct and which is uniformly bounded away from zero on Ed. (We may 
as well also assume that wt = 1 on S^.j Now, for any 5 G R, and any 4> £ C k ' a {^T), 
set 

\\<t>\\k,a,& = II^T 4>\\k,a, 
and letC$' a (T,T) denote the corresponding normed space. 



Proposition 7 Fix any 5 E R. For T sufficiently large, the mapping 

-,k- 

is an isomorphism. 



£ T :C 5 fc+2 '"(S T )^C^(S r ) 
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Proof: The action of Ct on C^ +2)OL (Y*t) is equivalent to the action of the conjugated 
operator w^ s Ctw t on the unweighted space C k+2 ' a . As 5 varies, these mappings are all 
Fredholm of index zero. Furthermore, any element of the nullspace of Ct is in any one 
of the spaces Cg +2 ' a . Therefore, to prove that any one of these maps is an isomorphism, 
it suffices to show that the nullspace of Ct is trivial when 5 = 0. In addition, by elliptic 
regularity, it even suffices to show that there is no element i] 6 C°(£t) such that Ct^j = 0. 

Assume this is not the case, so that there exists a sequence Tj — > oo and functions r/j £ 
C°(Y,Tj) satisfying CTjiVj) = 0- Write Ctj = Cj and normalize rjj so that sup \rjj\ = 1. As 
in the proof of Theorem ||, we shall show first that we may extract a nontrivial limiting 
function r], and then show that the existence of this function leads to a contradiction. 

Suppose first that the supremum of on T,* R remains bounded away from zero, 
at least for infinitely many j. Then using local elliptic theory to control higher order 
derivatives, we may extract a subsequence rjji which converges to some limiting function 
j] =k on £*. Since the coefficients of Cj converge on compact sets to those of C c , the 
linearization of the Lichnerowicz operator with respect to the initial data set (7 c ,/^c) at 
tfj = ij) c , we see that 

C cV = A lcV - i (r( 1c ) + 7|/ic| V c " 8 + y^c) V = 0. (22) 

The first two terms, A 7c — ^R(~f c ), constitute the conformal Laplacian for j c . Since 
7 C = V ; c~ 4 7) the familiar conformal covariance property of this operator leads to 

A 7c ?? - I fl(7c) V = € (a 7 (C^) " ~ ^(7) ■ (23) 

Using this in (|22"|), as well as the fact that |/i c | 7 = ip\ 2 \^ n and dividing through by ip^ 



we conclude that Cu = on £*, where u = ^/V'c- (Note that £ in this equation is the 

1/2 

operator which extends smoothly across the pj.) Next, r\ is bounded on S* and ip c ~ r- 

— 1/2 

near , and so, in Bj, \u\ < Cr j ' . This singularity is weaker than that of the Green's 
function in three dimensions, and so u extends to a weak solution of Cu = on all of 
S. But we have already remarked earlier that the operator C on £ has trivial nullspace, 
which means that u = 0. This is a contradiction. 

It remains to address the case in which the supremum of \rjj\ on Y,* R converges to 
zero. In fact, if sup \r]j\ is attained at some point qj € S-r,-, then we may apply the same 
argument as before if qj remains within a bounded distance of the end of the cylindrical 
piece CTj ■ If this distance tends to infinity, then choose a new linear coordinate s' which 
is centered at the point qj (so that s' = at that point) and is a translation of the 
coordinate s. Using local elliptic theory again we may extract a subsequence which 
converges on any compact set, along with all its derivatives, to some nontrivial function 
r] defined on the complete cylinder tx5 2 . Furthermore, the metric 7y. converges in this 
limit to the product metric 70 = (ds') 2 + h, R(jt-) ~^ 2, and the terms \/j,t -| 2 Vv 8 an d 
ip T , both converge to zero. Hence Cj converges locally on compact sets to A 7o — 1/4. 
Thus we have obtained in this limit a nontrivial bounded solution of the equation 

£ 7o r ? = A 7o r ? -^ = (24) 
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on the entire cylinder 1 x 5 12 , and such a solution obviously cannot exist. Hence this 
case too leads to a contradiction. 

This completes the proof of Proposition [?]. □ 

For T sufficiently large, let Qt denote the inverse of Ct acting on Cg' a . In other 
words, we have 

g T :Cj' a (E T )-^Cj +2 ' a (S T ), (25) 

with CtGt = Gt^t = I- Although of course Qt depends on the chosen weight 6, we 
suppress this dependence in the notation. 



Proposition 8 If < S < 1, then the norm of the operator Qt in (25) is uniformly 
bounded as T — > oo. 



Proof: The strategy is the same as before. If the result were false, then there would 
exist a sequence Tj — ► oo and functions fj G Cg' a (ETj) such that, writing Qj = Qt j7 

\\fj\\k,a,5 ->• and \\Qj(fj)\\k+2,a,S = 1 

as j — > oo. 

Setting Vj = Qj(fj), we argue exactly as in the preceding proof. Either some sub- 
sequence of the Vj converges to a bounded, nontrivial function v on T,* which satisfies 
C c v = 0, which leads to a contradiction, or else some subsequence converges on the 
cylinder. This second case is precisely why we have introduced the weight function w T .. 

Write wtj = Wj and Cy. = Cj. Since | \vj\ \k+2,a,6 = 1, we have 

1 > supsuptf^l^jl > c > 0. 

3 Cj 

Suppose that sup^. wj s \vj\ is attained at some point qj G Cj. We may also assume that 
\vj\ — * on any fixed compact set of X* away from Cj. Now, by assumption \vj\ < Wj, 
and if we renormalize by setting Vj = e Tj ^Vj, then this becomes 

\vj(s,9)\ < (cosh(s/2)) 5 

on Ctj, with equality attained at the point qj. Writing qj = (sj,9j), then we renormalize 
yet again, setting Vj(s',6) = Vj(s' + Sj,6)/(cosh(sj/2)) s . Hence 

V cosh 2 S j ) 

with equality at the point (s',9) = (0,9j). In particular, Vj(0,9j) = 1. Now we can pass 
to a limit. If Sj — * ±oo, then the right hand side here converges to e 5s I 2 or e~ 5s / 2 , while 
if Sj remains bounded, then the right side converges to a function which is bounded by 
C(cosh(V/2))' 5 for some C > 0, while the left side converges to a nontrivial function v 
which satisfies C,~ v = 0. But any solution of this equation blows up at least like e' s 'Z 2 
either as s' — > oo or as s' — > — oo. Since < S < 1, this rate of growth is incompatible 
with the previous inequalities, and so we reach a contradiction once again. □ 
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6 Proof of the main Theorem 

We wish to solve 

N t {^t + vt) = 0. (26) 
Expanding Mt{^t + ?]) in a Taylor series about r\ = we have 

M T {rh + v) = E T + Cr(ri) + Qt(v), 
where Er and Ct are given by (|i~8|) and (|2l|) respectively, and 



Qt(v) = Mrifh + v) -E T -£ T (v) 

= \\M 2 1T {{i> T + v)~ 7 -i>T 7 + 7 ^ s v) 

-^r 2 (lOV'rr? 2 + W^W + + V 5 ) ■ 

is the "quadratically vanishing" nonlinearity. If we require that 77 G C^ +2 ' a (T>t) satisfy 
|?7| < c0t for some constant c < 1, then from the expansion above we easily obtain an 
estimate of the form 

llQTmik^a^CMt^ (27) 

for some constant C independent of T. Equation (^6|) is equivalent to 

VT = ~Qt (Et + Qt(vt)) 

where Gt is the inverse of Ct acting on C^ +2 ' q (St). We will find a solution to this 
equation in a small ball about the origin in C^ +2,a (^x) for $ G (0)1)- To determine 
the appropriate size of the ball in which to work, note that from Proposition |(| for any 
5 G (0, 1), we have 

\\E T \\k,a,S < C e" T / 4 (28) 
for a constant C which is independent of T. This suggests the following 
Definition 4 Fix k G N, a G (0, 1) and 5 G (0, 1), and let v G M + . T/ien we define 

= {ti G c£ +2 ' a (£ T ) : IMIfc+W < ^e" T / 4 }. 
Thus, for any 77 G M u , we have 

M < ve~ T ^ on Sj- \ Ct and |t/| < v e~ ST ^ipT on Ct- 

In particular, when T is sufficiently large, the quadratic estimate (^7j) is valid in M u . 
Using Proposition || together with ( p8| ) we have 

||g T (^T)IU + 2, a ,5<CMe- T / 4 (29) 

where M is the bound on Qt acting on C^'"(St). In light of this we choose v = 2CM so 
that we have Gt(Et) C ^ u /2- The following Lemma follows easily from the fact that we 
have the quadratic estimate above on Qt(v) f° r 1 SB„. 
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Lemma 1 For v = 2CM as above and T sufficiently large we have for any 771 , 772 £ B„ 

||2tM - QT(m)\\k+2,a,S < T^jWVl - m\\k+2,a,5- (30) 
We can now establish the following 

Proposition 9 For v as above and T sufficiently large the map 

X] .— > T( v ) = -G T (E T + Q T ( V )) 
is a contraction mapping on M u C C^ +2,a (T,t) ■ 

Proof: Let 77 £ B v be given. Taking 771 = 77 and 772 = in (^) we see that —Qt(Qt(v)) e 
B„/2. This, together with our choice of ^ shows that T : B y — ► B„ as required. Moreover 
for 771 , 772 £ By we have 

||T(tti) - r(77 2 )|| fe+2)Q)l5 = \\g T (Et + QHm)) - Qt (e t + Q T (m)) lU+w 

< m|| Q T (m) - Qr(m) II 

1„ 

< ~ m\\k+2,a,S- 

This completes the proof of Proposition ||. □ 
Theorem [l] from § L2 is now a direct consequence of the following. 



Theorem 3 The Lichnerowicz equation (26) has a unique solution 4>T = ipT + VT with 

Proof: This follows immediately from the contraction mapping theorem, together with 
Proposition 0. □ 



7 Modifications when E is complete but noncompact 

Thus far we have assumed that the three-manifold S is closed. However, it is also 
interesting from a physical perspectical to consider gluing solutions of the constraint 
equations on various classes of complete manifolds. Of particular interest are the two 
cases when (£,7) is either asymptotically Euclidean or asymptotically hyperbolic. For 
each of these cases, there are well-established results concerning the existence of CMC 
solutions of the Einstein constraint equation ]7j , [ 14 1 , [Q] , [|| . In this section we show that 



with only minor modifications, the gluing constructions developed here may be extended 
to these other settings. Since most of the modifications are common to both cases, we 
begin with a more general discussion which describes the alterations and new ingredients 
needed to extend the arguments, modulo a few specific analytic facts. We follow this 
with precise statements of definitions and facts needed for these two particular cases, 
along with appropriate references. 



27 



7.1 General preamble 

Looking back, there are only a few places in the constructions and arguments of the 
preceding sections where truly global issues arise; elsewhere the argument localizes in 
the neck region. Specifically, the global structure of £ enters the analysis only when 
establishing the existence of the inverses Gt and Qt for the vector Laplacian Lt and 
the linearized Lichnerowicz operator Ct, and also, more mildly, when determining the 
dependence on T of the norms of these inverses. 

Thus suppose that (£,7) is a complete, noncompact, three-manifold, and (7, II) is 
an initial data set on £ with trLT = r constant. In the noncompact setting we do not 
require that LT ^ 0. Fix points pi,P2 £ £ an d balls Bj 3 pj and construct the approximate 
solutions (7TiA*T)V ; t) as in §2. Then proceed to the main part of the argument, where 
first [it is modified to be transverse-traceless with respect to 7t using the operator Lt, 
and then tpT is modified so as to be a solution of the Lichnerowicz equation. 

The function spaces we shall use in the two cases of interest are slightly different, 
and so we phrase things more abstractly in this preamble. Thus suppose that there are 
Banach spaces E and F of vector fields and E' and F' of functions - these will be familiar 
weighted Holder spaces in practice - such that the mappings 

L : E — ► F, C:E' — ► F' (31) 

for the vector Laplacian L and the linearized Lichnerowicz operator C on (£, 7) are both 
isomorphisms. 

Remark 3 In the compact case, we assumed that there were no conformal Killing fields 
vanishing at both p\ and P2, and this was enough to give the invertibility and bounds for 
the inverse as T — > 00 for Lt; likewise, the assumption that II ^ was used to prove 
invertibility of Ct- In both the asymptotically Euclidean and asymptotically hyperbolic 
settings, neither assumption is required because we shall use spaces of functions and vector 
fields which decay at infinity, and for these it is not difficult to obtain invertibility for Lt 



and Ct directly. Thus the restriction that the maps (31) are both isomorphisms is not 
too limiting. 

We must also assume that the vector fields X £ E decay sufficiently so that the 
integration by parts 

(LX,X) = (VX,VX) (32) 

is valid for all X € E. 

Since we shall always use function spaces which are localizable, i.e. preserved by 
multiplication by functions in C(J° , we may define new function spaces on £y by patching 
together vector fields or functions in any of these spaces with vector fields or functions 
in the appropriate Holder spaces on Ct (including functions measured with the weight 
factor w s as required in Propositions 7 and 8). We label these spaces E(Y±t), F(T,t), 
E' s (Y<t) and Fl(Ej>). It is clear that, presuming as above that L and C are isomorphisms, 
then 

Lt : E(Er) — F(Er), C t : E' S (Z T ) — i^(Sr) 
are Fredholm for any T. For simplicity, we shall usually drop the £t from this notation. 
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Proposition 10 The operator Lt '■ E — ► F is invertible when T is sufficiently large. 
The norm of its inverse Gt, as a mapping F — > E, is bounded by CT 3 as T — > oo. 



Proposition 11 When 5 £ (—1, 1), the operator Ct '■ Ei — ► F$ is invertible, with 
inverse Gt uniformly bounded as a mapping Ft —* E' s as T —* oo. 

We shall prove only the first of these two results; the second is proved in exactly the 
same way, but in fact is even easier to verify. 



Proof of Proposition 10: Let us suppose the result is false, so there is a sequence of 



values Tj — > oo, and vectors Xj and Yj such that Lt Xj = Yj, but with = 1 and 

II^'Hf = o(T -3 ). Also, write E(Y*Tj) and F(Y,Tj) simply as E and F. 

Following the argument of Theorem 0, let qj denote a point where ||Xj||_e attains 
its maximum. There are three cases to consider: the first is when qj lies in Cy. and 
the norm of Xj on any compact set of £* tends to zero, the second is when qj remains 
in some fixed compact set of £*, and the third is when qj leaves every compact set of 
£. In the first two cases, the argument proceeds exactly as in §3.2 and §3.3; either we 
produce in the limit an 'illegal' solution X on the complete cylinder, or else we produce a 
nontrivial vector field X on £* such that L C X = there. Here L c is the vector Laplacian 
for the metric 7 C on X* which agrees with 7 outside of the balls Bi, £ = 1,2, and has 
asymptotically cylindrical ends near these points. This solution X is bounded along these 
cylindrical ends and it lies in the space E on X* \ (B\ U B2). In order to carry out the 
usual integration by parts, to conclude that T> lc X = 0, we use the analysis from §3.2 and 
the assumption (^) for X £ E. Hence X is conformal Killing, first with respect to 7 C , 
and hence with respect to 7, on X*. As in §3.2, it extends to a smooth solution across 
the marked points pi, and therefore gives a nontrivial element X £ E of the nullspace of 
L on all of X, which is a contradiction. Thus it suffices to handle the third case. 

Suppose then that \Xj\ decays to zero on every compact set of X. Let x be a smooth 
cutoff function which vanishes near the points pe and equals one outside the balls Bg. 
Then Xj = \Xj satisfies L^Xj = xYj + [Lt^, x\Xj = Zj. Both terms on the right decay 
to zero in the space F on X (not just X*), whereas clearly ||Aj||£ > c > 0. However, these 
vector fields are all defined on (£,7), and GZj = Xj. Since G is a bounded operator, 
this is again a contradiction. □ 



7.2 Asymptotically Euclidean initial data 

We now specialize the discussion above to the case in which (£,7) is asymptotically 
Euclidean and (7, n) is an appropriate data set on this manifold. 

Recall that a metric 7 is said to be asymptotically Euclidean (AE for short) if each 
end of X is diffeomorphic to the exterior of a ball in M 3 , and in the induced coordinates 
the metric 7 decays to the Euclidean metric at some rate. More precisely, we require first 
that 

$:£\iT«|J(R 3 \^.), 

where the disjoint union on the right is finite. For simplicity we assume that £ has only 
one end. Then, we assume that in the naturally induced Euclidean coordinates z, 
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along with appropriate decay of the derivatives. 

To formalize this, suppose that E \ K as R 3 \ Br q , and then define 



Aae(S) = {u : SUp|li| < OO, | \u\ \o,a,K < OO, 



\u(z)-u(z')\R a 
sup sup j < OO 

R>R z^z> \Z-z'\ a 

R.<\z\,\z'\<2R 



Ai' B Q (S) = {n:||n|| fcia ^<oo, (1 + \z\ 2 )^ dgu G A°£, for \(3\ < k}, 



r U A^(T,) = {u G A A '° : u = r u v on E \ K and u G A^'" there} 

Thus, roughly speaking, a function u is in r^A^ provided it is in the ordinary Holder 
space A k,a on any fixed compact set in E, and if \dzu\ < C^|z| _1/_ '^' along the end(s), 
for < k, with an appropriate condition for the Holder derivatives of order k + a. 
We now define 

M k _!^(T,) = {metrics 7 : 7^ - <% G r~ u Aj£ }. 

(We are abusing notation slightly here since the metric coefficients 7^ are only defined 
on the end E \ K.) Such metrics have been much studied in relativity, and Q collects 
a number of analytic results which are appropriate for the study of the Laplacians and 
other natural geometric operators associated to these metrics. 

In this AE setting it is natural to focus on initial data sets which would evolve into 
asymptotically Minkowski spacetimes, and so we assume that r = 0. In addition, since 

7 G M fc + 2 ' a (E) =^ Riem 7 G r-^A^E), 

the Hamiltonian constraint equation (Q) suggests that we assume that 

/x G r-l-^E). 

There are many other possible choices of function spaces in which to work. The most 
customary ones in the literature surrounding this problem are weighted Sobolev spaces. 
We could certainly have used these here too, but because weighted Holder spaces have 
been used everywhere else in this paper it seems more natural to use spaces of this type 
here as well. 

We now define the class of elliptic AE operators. Suppose that P is an elliptic (second 
order) operator on E; we shall assume that there is a constant coefficient elliptic operator 
P of order 2 on E \ K ra M 3 \ B Ro such that 

P = P 0+ ^2 ap(z)dg, where a p G r~ v A k A t 2 ~ m ' a . 

I/3K2 
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Thus by construction, if P is of this type, then 

P : r- u A k A t 2 ' a (Z) — ► r- u - 2 A k ^(Z) (33) 

is bounded. There are additional conditions required to ensure that this mapping is 
Fredholm. These are determined by the indicial roots: the number A is an indicial root 
of P (on one end of E) if there exists a function 4>(0), 6 = z/\z\, such that P(\z\ X (j)(9)) = 
0{\z\^~ 2 ~ v ). Notice that no matter the value of A and choice of <fi, the right side is 
always of the form Pq{t~ v <j){6)) + 0(\z\^~ 2 ~ u ), where the first term here is homogeneous 
of degree A — 2. Hence indicial roots are determined solely by the 'principal part' Po of 
P. Continuing, if we express r 2 Po i n polar coordinates as 

i+|/3'|<2 

(where the coefficients aj^i are constant) then A is an indicial root for Po an d hence for 
P if and only if 

£ a^,A^>(0) = O 

j+l/3'l<2 

for some function (j)(9). In particular, there is a denumerable set of indicial roots {Aj}, 
and these are related to eigenvalues of the 'angular part' of the operator Pq. 



It is known [23], y] that if — v ^ {Re Aj}, then the mapping (33) is Fredholm (whereas 
if — v is the real part of an indicial root, then this mapping does not have closed range). 
We are concerned with the question of whether this mapping is an isomorphism, or at 
least surjective, and if so, for which range of weights, for the particular operators L and 
C. It is clear from the asymptotic flatness of the data that the principal parts of L and 
C are 

11 

Lq = — A H — Vdiv and Cq = A, 
2 6 

respectively. (The Laplacian on the right in this expression for Lq is the Laplacian on 
vector fields on R 3 , which is the same as the scalar Laplacian applied componentwise 
with respect to the standard Cartesian decomposition.) From here we compute that 

• the set of indicial roots for both L and C is Z. 

For the scalar Laplacian this is well known and derives ultimately from the fact that the 
global temperate solutions of Cqu = are polynomials, and these have integer rates of 
growth. The negative indicial roots arise from considerations of duality between weighted 
L 2 spaces. Exactly the same reasoning may be applied to determine the indicial roots of 
Lq, for it is easy to see using the Fourier transform that the only global temperate vector 
fields which satisfy LqX = on R 3 must have polynomial coefficients. 
Then we have the following 



Proposition 12 When P is equal to either L or C, the mapping (35) is Fredholm of 
index zero when — 1 < — v < 0. 

We remark that the fact that these mappings are Fredholm follows from the results 
quoted above and the fact that when — v is in this range it is not an indicial root of 
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either operator; the fact that they have index zero can be derived from the symmetry 
of these operators. When —v is not in this range, and is not an indicial root, then 
the mapping ( |33"| ) is still Predholm, but no longer of index zero, hence cannot be an 
isomorphism. 



Proposition 13 For P = L or £ and —1 < —v < 0, the mappings in (35) are isomor- 
phisms. 

Proof: It suffices to show that neither L nor C have elements in their nullspace which 
decay like r~ v as r — ► oo. For the linearized Lichnerowicz operator C this follows directly 
from the maximum principle. For the vector Laplacian L, we first note that if X S 
r~ v A 2 ^ (£) with —1 < —v < 0, then the boundary term in the integration by parts 
(LX,X) = ||2?X|| 2 vanishes, and so T>X = 0. We then use a result of Christodoulou 
§ which says that there are no conformal Killing vector fields on AE manifolds which 
decay at infinity. (In fact, Christodoulou's theorem is stated only for AE manifolds which 
are globally diffeomorphic to M 3 , but one readily verifies from his proof that the result 
extends to general AE manifolds). □ 
In summary, we have proved 

Theorem 4 Let (£,7,11) be an asymptotically Euclidean initial data set, as described 
in this section; in particular, r = ir 7 (II) =0 on £. Choose any two points p\,P2 & £, 
and suppose that if either of these points lies on a compact (closed) component o/£, then 
there does not exist a conformal Killing field X on that component which vanishes at that 
point, and furthermore, II ^ on that component. Then there exists a one-parameter 
family of asymptotically Euclidean CMC initial data sets (Tt, Ht) on £t satisfying the 
same estimates as in Theorem 0. 



7.3 Asymptotically hyperboloidal initial data 

Let us now turn to the case where (£,7) is asymptotically hyperbolic (or hyperboloidal, 
which is the terminology favored in the physics literature). We work within the framework 
of conformally compact manifolds, which is quite natural from Penrose's formulation. 
This is explained carefully in the introduction to Q. 

Recall that a metric 7 on the interior of a manifold with boundary £ is said to be 
conformally compact if 7 = p™ 2 7, where 7 is a nondegenerate metric on the closure £ 
and p is a defining function for <9£, i.e. p is nonnegative and vanishes precisely on the 
boundary, with nonvanishing differential. We shall assume that 7 and p are both smooth, 
or more accurately, polyhomogeneous (i.e. having complete expansions involving powers 
of a smooth defining function as well as powers of its logarithm); this latter assumption 
is most reasonable because it occurs generically |^]. Though we shall not do so, it is 
also possible to assume that 7 and p are in C fc+2 ' Q (£), and only minor and obvious 
perturbation arguments are needed to handle this more general case. 

The metric 7 is complete on the interior £ of £; moreover, assuming that \dp\~ = 1 
at p = 0, then the curvature tensor Riem g is asymptotically isotropic as p — ► 0, with 
sectional curvatures tending to —1. We call any such space AH (standing for either 
asymptotically hyperbolic or asymptotically hyperboloidal, depending on the reader's 
preference). As for the tensor II = p + 5T7, we assume that r is constant, of course, and 
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also that the tensor coefficients = 0(p~ 1 ), or equivalently // J = 0(p 3 ). (In either 
case, we are writing the components of the tensor [i relative to a smooth background 
coordinate system on £.) Assuming that r = 3 and with the curvature normalization 
above, the constraint equations become 

div 7 (/i) = 0, i? + 6 = |/i| 2 . 

(Recall that |/i| 2 = 0{p 2 ).) 

The existence of AH initial data sets (7, II) satisfying these hypotheses, as well as the 
polyhomogeneous regularity of this data (assuming that 7 is smooth) is the main topic of 
Q. The special case of these results, in which one assumes that II is pure trace, is already 
handled in H. Some of the main tools required here are based on the linear analysis of 



geometric elliptic operators on conformally compact spaces as developed in [24|, although 
the authors of [|j and Q base their work on the somewhat later and less precise, although 
simpler methods of cf. also [22| for a sharper version of these methods. The regularity 



theory of U is very closely related to the results of |2J 

Quoting the main result of Q , we assume that £ is endowed with an AH initial data 
set (7, n), where 7 and \i are polyhomogeneous. That is, with 7 a smooth metric on £ 
and x a smooth defining function for <9£, we have 



and 



fc>0 i=0 

fc>o e=o 



x 



where all coefficients a^, fi[j , are smooth on S and for each j, vanish when 

k > Nj. (In addition, a 2 e = for £ > 0.) 

We require the mapping properties of the vector Laplacian L = — div o T> and the 
linearized Lichnerowicz operator C = A — \p\ 2 — 3 associated to this data. These follow 



directly from the general theory of p4[ , cf. also [22]. 

Choosing local coordinates (x,y) near <9£, then recall from these papers that any 
operator of the form 

P= a jf3(x,y)(xd x y(xd y ) a 

j+M<2 

with coefficients smooth (or just polyhomogeneous) is said to be uniformly degenerate 
(of order 2). Its 'uniformly degenerate symbol' is defined to be 

i+l/3|=2 

and P is elliptic in this calculus if this symbol is nonvanishing (or an invertible matrix, 
if P is a system) when (£,77) 7^ 0. It is not hard to check that both of the operators L 
and C above are elliptic uniformly degenerate operators. 

The function spaces on which we let these operators act are the weighted Holder 
spaces x A££(E). These are taken to be the usual Holder spaces on a compact subset of 
£ and, in a boundary coordinate chart are defined as follows: 
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A AH ( s ) = i u ■ SU P M < °°> SU P 71 jr—\ < 00 

1 (x-x' + y-y \r 



A*£(E) = {u : (xd x y{xd y fu E A°£, j + < fc}, 



x"A a 'h (S) = {u = x v v : v E A a 'h }• 



Clearly 



P : x-A A + 2 '"(S) — > x^A^(S) 

is bounded for any v, but just as in the AE case, restrictions on v must be imposed to 
ensure that this map is Fredholm. 

To state these, we say that a number A is an indicial root for P if Px x = C(x A+1 ) 
as x — > 0. It is always true that Px x = 0(x x ) and the indicial root condition requires 
that A satisfy a polynomial equation involving the coefficients Ojo(0,y). (It is possible 
that these roots could depend on y, but in our application this does not occur.) A brief 
calculation, using that \n\ 2 -> as x — > 0, gives: 

• The indicial roots for L are and 4, 

• The indicial roots for C are —1 and 3. 



We now state one of the main theorems from [24|, specialized to this particular setting: 



Proposition 14 Suppose that < u < 4 and — 1 < 1/ < 3. Then the mappings 

L : x^+ 2 < a (£) x"A^(£) 

and 

£ : x^'A^X) — > x^'A^(S) 

are Fredholm of index 0. Furthermore, the nullspace of either of these operators is inde- 
pendent of v E (0,4) and z/ E (—1,3). 

To conclude that these maps are isomorphisms, it suffices to prove 

Proposition 15 For v and v' in the ranges stated in the previous proposition, the maps 
L and C are infective. 

Proof: For C this is quite easy. In fact, if L<\> = 0, then < Cx v for any u' E (—1,3), 
and therefore (j) — > at 9S. But the terms of order in C are strictly negative, and so 
the maximum principle implies that <f> = 0. As for L, we note that if LX = 0, then the 
coefficients of X are 0(x 4 ), and so we may perform the usual integration by parts to 
obtain that T>X = 0. By the conformal invariance of the nullspace of T> = 2? 7 , we also 
have that T>—X = 0. Thus X is a (smooth or polyhomogeneous) conformal Killing vector 
field for the metric 7 on S, with coefficients vanishing to order 4 at the boundary. Since 
we already know that X is polyhomogeneous, the equation T>— X = now implies that 
X vanishes to infinite order as x — > 0, and it is known that no such conformal Killing 
field can exist [Hi- Thus we have shown that 
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Proposition 16 The mappings in Proposition 1^ are both isomorphisms. 



Therefore, we may use the spaces occuring here as the spaces E, F, E' and F' needed 
in Propositions 10 and This proves 



Theorem 5 Suppose that (S, 7,11) is an asymptotically hyperbolic initial data set, so 
that j is a conformally compact metric and r = £r 7 (II) = 3. Let pi,P2 £ S, and suppose 
that if either of these points lies on a compact (closed) component o/S, then there does 
not exist a conformal Killing field X on that component which vanishes at that point. 
Then there exists a one-parameter family of asymptotically hyperbolic CMC initial data 
sets (Pr,IIx) on the manifold T>t satisfying the same estimates as in Theorem^. 



8 The geometry of the neck 

We have now shown in the three main settings (S compact or AE or AH) that it is 
possible to produce one-parameter families of CMC initial data sets (Tt,Ht) which are 
arbitrarily good approximations to the original initial data set (7, II) away from the points 
pi, p2- As T — > 00, various geometric quantities associated to these tensors degenerate 
in the neck region. We describe this in more detail now. 

We shall only consider quantities which depend on at most two derivatives of the 
metric, and so we fix any k > 2 and a £ (0, 1) and construct the initial data sets as in 
§§2-7. 

Recall the notation = £ \ (B R (p{) U B R (p 2 )). Next, for b £ (0,R), let Aj(b) = 
Br{Pj) \ Bb(Pj), and finally denote by Cx,b the part of the neck region not covered by 
these annuli. Then for each T > 0, 

S T = U (Ax(b) U A 2 {b)) U C T , b - 

In either of the annuli and in the truncated neck region we shall use the (s, 9) coordinates, 
so that, for example, A^b) = {T/2 - (5 < s < T/2} x S 2 and C Ti& = {-T/2 + /3 < s < 
T/2 — /?} x S 2 , for = — log b. Write tp^'jT = It] of course, this metric is the same as 7 in 
T,* R , and the deviation between these two metrics in either of the half-necks < \s\ < T/2 
is 0((e~ T cosh s) 2 ). We shall estimate \Tt~ 1t\^ t (^- e - the difference between the actual 
solution IV and the approximate solution jt) in each of these regions. 
By definition, Tt = (tpr + ??t) 4 7t- Hence 

I It - 7r| 7T ~ ^ 4 ((V"t + ??t) 4 - tpr) < C\vt/iPt\- 

Recall also from §|6] that \rjr\ < ve~ T l^w & Tl Furthermore, 1/2 < wx/ipT < 2 everywhere. 
Therefore: 

• In ~E* R , 7t = 7 and ipT ~ 1- Thus 

|r T - 7t|7t < Ce- T ' 4 (34) 

in this region. 

• In all of Ct, we have 

|r T - 7t| 7t < Ce- T / 4 (e- T / 4 cosh(s/2)) <5 - 1 = C e - <5T / 4 (cosh(s/2)) 5 - 1 . (35) 
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• In particular, in Ai(b), \rfr\ < ve T//4 ui T , and so 

|r T - 7t|t t < Ce-^e-^V/ 2 ) 5 " 1 = Ce-^e^ 1 ^ 2 < Ce^W 2 ^/ 4 . (36) 
The estimate in A2(b) is identical. 

Notice in each of these cases that since 5 < 1, this difference of tensors is much smaller 
than 7t, and therefore geometric quantities relative to Tt are very well approximated by 
the same quantities relative to jt- 

The following two results are immediate corollaries. 

Proposition 17 Let c be any C 1 path in £t which runs through the neck region. For 
simplicity we shall assume that it has the simple form s i— > (s,9q), —T/2 < s < T/2. 
Then 

length rr (c) ~ length^ T (c) + 0(e" 5T / 4 ). 

In particular, for the portion of c which lies in either of the annuli Aj (b) , we may sub- 
stitute length 7 (c) for lengthy (c) on the right side here. 

Proposition 18 In CT,b, the full curvature tensor Riemr T satisfies 

|Riemr T |7 T < < Ce T (cosh s) -2 . 

In particular, the scalar curvature Rr T satisfies this same estimate. 

Of course, we could also compute finer asymptotics for the full curvature tensor and for 
the scalar curvature function. 

Estimates for the extrinsic curvature tensor are also easy to obtain. We have already 
shown that the deviation of the trace-free part of TIt from that of II on T,* R is bounded 
by T 3 e -3T//2 . On the other hand, relying on the fact that |/Li| 7c < Ce~ T / 2 coshs together 
with the preceding estimates, we see that 

|n r |r T < Ce T/2 (coshsr\ on C T . 

We may use this as a check on this last proposition: the constraint equation (||) gives 
the scalar curvature in terms of |I1t| 2 t , and we see that the estimates are the same. 
However, we also see from this that if the original tensor \i = in both balls B\ and B2, 
then |ny|r T is bounded in Ct, and hence Rr T is also. 

9 Examples and applications 

There are a number of constructions and applications of special interest in general rel- 
ativity to which our results pertain. In this section we briefly describe several of these. 
This discussion is intended to be informal and descriptive, and so we do not present these 
results formally, but given the rest of the contents of this paper, it should be obvious how 
to do so. These examples are separated into different subsections, but these divisions are 
artificial and are intended for guidance only. 
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Adding wormholes 

The first construction is that of adding a physical wormhole to (almost) any globally 
hyperbolic solution of the vacuum (Lorentzian) Einstein equations which admits a CMC 
Cauchy surface. One proceeds by: a) choosing such a Cauchy surface S; b) choosing any 
pair of points on S; c) verifying the conditions of Theorems [l], ^| or ||; d) carrying out 
the gluing construction as outlined in §§2-6; and finally, e) evolving the resulting initial 



data. Standard well-posedness theorems |ll|], [|13[] guarantee that there is a unique (up 



to spacetime diffeomorphism) maximal spacetime development of this initial data set. In 
other words, there is a spacetime which satisfies the vacuum Einstein equations, has the 
topology Tit x restricts to the wormhole initial data on Et x {0}, and also contains 
every other spacetime which satisfies these three properties. Of course, this spacetime will 
probably not last very long, and indeed it is likely that the evolving wormhole will pinch 
off quickly, thus preventing the spacetime from evolving further in time. Nonetheless, 
there is a spacetime solution which evolves from this initial data and which therefore 
contains a spacetime wormhole connecting the chosen points; in addition, the evolving 
spacetime remains essentially unchanged from its previous unglued state outside the 
domain of influence of the initial data on the wormhole region. 

There do exist spacetimes for which this wormhole addition cannot be carried out: 
these are the ones which do not carry any set of CMC data which matches the hypothe- 
ses of our theorems. A simple example of this phenomenon is the 'spatially compact 
Minkowski space' T 3 x R, since all CMC slices of this spacetime are compact and have 
vanishing extrinsic curvature, and hence the hypotheses of Theorem 1 are not satisfied. 

Multiple black hole spacetimes 

Over 40 years ago Misner |3(| constructed a number of explicit initial data sets in- 
volving two asymptotically Euclidean regions connected by a wormhole. Soon thereafter, 
he and others, cf. [17], recognized that since the wormhole inevitably contains a mini- 



mal two-sphere, and since the presence of such a two-sphere necessitates the existence of 
an apparent horizon, the spacetime developments of these initial data sets must contain 
black holes. In subsequent work, Misner constructed explicit (series form) sets of 
initial data with more than one wormhole; the spacetime developments of these data 
must then contain more than one black hole. Misner's data has been used extensively by 
numerical relativists (TTJ] in studying multi-black hole spacetimes. 

Our gluing constructions yield a much wider class of multi-black hole solutions. Most 
importantly, we can successively add black holes to any given asymptotically flat space- 
time. (We say that a spacetime is asymptotically flat if it is foliated by asymptotically 



Euclidean initial data sets; see chapter 11 in [34| for a more rigorous definition.) The idea 



is as follows: start with any asymptotically flat spacetime which admits a maximal (i.e. 
zero mean curvature), asymptotically Euclidean hypersurface. Choose such a maximal 
Cauchy surface with AE data (£,7, II) and a point on this hypersurface, and then glue 
on a copy of Euclidean space M 3 with zero extrinsic curvature (which is a maximal data 
set for Minkowski space). The result is a new initial data set with two AE regions con- 
nected by a wormhole, with the data largely unchanged on either side of the wormhole. 
As before, the wormhole contains a minimal two-sphere, and so black hole formation is 
inevitable in the spacetime development of this data. We may now continue this process, 
taking this new data and gluing on a further copy of maximal data for Minkowski space. 
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This yields a spacetime with two black holes. Proceeding further, we obtain spacetimes 
with an arbitrary (finite, or even infinite) number of black holes, with each additional 
black hole leaving the data for the previous black holes largely unchanged. 

If one glues two AE sets of data together, the resulting data set has two AE ends. From 
a mathematical point of view, there is nothing to distinguish one of these ends from the 
other. However, in using Einstein solutions to model spacetime physics, one often makes 
a choice of one of the ends, and its spacetime development, to be the physically accessible 
"exterior" , with the other end corresponding to a physically inaccessible "interior" . The 
exterior is the home of the "observers at infinity" while the interior is somewhere "inside 
the black hole" . Apart from considerations of physical modelling, the choice is arbitrary. 

Multi-black hole spacetimes have several ends, and any one of these may be regarded 
as the exterior. However, to model spacetimes with observable interacting black holes, 
one must construct multi-black hole initial data sets with the chosen exterior asymptotic 
region directly connected to each of the wormholes. Each of these wormholes leads 
through to its opposite interior end. These interiors may be independent of one another, 
or they may coincide as a single asymptotic region. In other words, one may envision the 
manifold as having, for example, two asymptotic regions connected by several necks. This 
case corresponds to Misner's "matched throat" data sets and is no more difficult to 
construct from our point of view than a solution with independent interior ends. 

We emphasize that the multi-black hole data sets we are describing here are not 
presented as explicitly as Misner's examples, because to obtain them one must solve the 
constraint equations in conformal form. In practice, there are efficient ways to do this 
numerically. As already noted, one virtue of our results is that we obtain solutions which 
leave the spacetime essentially unchanged away from the domain of influence of the added 
black hole. This feature is quite interesting and should be useful to numerical relativists. 

Adding AE ends to closed initial data sets 

Every closed three manifold £ admits a CMC solution of the vacuum constraint 
equations: we choose 7 to be a metric of constant negative scalar curvature #(7) = —6 
(these always exist), and set n = 7; (S, 7, 7) is then a solution of the constraint equations, 
with constant mean curvature r = 3. It is of interest to know whether every manifold 
of the form £ \ {p} with £ closed admits AE initial data sets. For a manifold £ which 
admits a metric of positive scalar curvature with no conformal Killing fields and with a 
nonvanishing symmetric transverse-traceless 2-tensor, our methods show that this is the 
case. To see this, we first use [18] to get that any such £ admits a solution (£,7, n) of 



the constraint equations with r = 0. Then by Theorem [|, we may glue the standard AE 
initial data set on R 3 with vanishing extrinsic curvature to (£,7, II) to get an AE initial 
data set on £ \ {p}. Note the closely related result, that when £ is compact, £ \ {p} 
admits a complete scalar flat metric, i.e. a time-symmetric AE initial data set, if and 



only if £ admits a metric of positive scalar curvature, cf. [25]. 

On the other hand, suppose that £ has the property that any metric of nonnegative 
scalar curvature on it is flat (for example, £ = T 3 has this property). Then by @, any 
initial data set (£,7, n) with r = has nonnegative scalar curvature. By the hypothesis 
on £, this means that 7 is flat, and so, using (§) again, we must have n = 0. But (£, 7, 0) 
does not satisfy the hypotheses of Theorem ^, and so we cannot use our gluing theorem 
to produce AE initial data sets on £ \ {p}. 
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We now comment briefly on the physical interpretation of the solutions we obtain 
by gluing AE ends to closed initial data sets as above. Suppose we view the compact 
side as the exterior, accessible region. This data has a minimal two-sphere as before, and 
therefore most likely evolves into a spacetime with an apparent horizon. But the problem 
is that there is no clear and unambiguous definition of a black hole relative to an exterior 
spacetime region which is not asymptotically flat. This makes the interpretation of the 
spacetime development in this case less transparent, but it also makes such spacetimes 
rather interesting. Indeed, one might be able to learn a lot about how to work with black 
holes in spatially closed spacetimes by careful study of these examples. 

As noted above, if a spacetime contains a wormhole connecting a pair of regions, one 
may choose either of the pair of regions as the exterior accessible region. The choice is 
mathematically arbitrary, but physically important for modelling. Thus the spacetimes 
one obtains by gluing together data from an asymptotically flat spacetime and data from 
a spatially compact spacetime, and then evolving, may be viewed either as a spatially 
compact cosmology with a black hole-like object in it (with the closed region as the 
exterior), or as a black hole with a nontrivial interior solution. The construction is the 
same; all we are doing is redefining which region we consider to be the exterior, physical 
region. The wormhole may cut off very quickly and render the interior solution physically 
irrelevant; this is the nature of black hole physics, and we believe that it can be studied 
meaningfully with these gluing constructions. 

Adding AH ends to closed initial data sets 

In distinction to the AE case, since AH initial data sets have nonvanishing mean 
curvature r = 3, it is possible to glue an AH initial data set to an appropriately chosen 
intial data set on any compact manifold E. In fact, as above, let 7 be a metric of constant 
negative scalar curvature R = —6 on E, so that (E, 7, 7) is an initial data set also with 
r = 3. By Theorem || we may glue to it any AH initial data set (in particular, the 
hyperboloid (H 3 , 70, 70), where 70 is the standard hyperbolic metric) at any point p G E. 
This proves in particular that for every closed manifold E, the punctured manifold E\{p} 
admits an AH initial data set. In fact, we may iterate this procedure and obtain AH 
initial data sets on E \ {pi, . . . ,pk} for any collection of distinct points on any closed 
3-manifold E. We note again that the analogue of this question for AE vacuum initial 
data sets is still open. 

As a simple and interesting example of this last construction, let (E,7) be a closed 
hyperbolic 3-manifold. By adjoining to it the hyperboloid as above we obtain an AH 
initial data set on E\{j?}. In fact, we may find such initial data sets which are arbitrarily 
small perturbations of the initial hyperbolic metric away from the neck. 

It is also straightforward to glue together two or more AH data sets. Not surprisingly, 
it is impossible to glue CMC AH data sets to CMC AE data sets because the constant 
mean curvatures in the two cases are different. 

We are left with the question of whether the spacetimes which result from evolving 
these glued AH initial data sets develop black holes. Again, if we view the spatially 
compact part of the data as the exterior, physically accessible, part of the solution, then 
this question does not have a definite meaning, since as before, black holes are only 
well-defined in spacetimes with asymptotically flat ends. 
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